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Abstract 

The sub-cycle dynamics of electrons driven by strong laser fields is central to the emerging field of attosecond science. We 
demonstrate how the dynamics can be probed through high-order harmonic generation, where different trajectories leading to the 
same harmonic order are initiated at different times, thereby probing different field strengths. We find large differences between 
the trajectories with respect to both their sensitivity to driving field ellipticity and resonant enhancement. To accurately describe 
the ellipticity dependence of the long trajectory harmonics we must include a sub-cycle change of the initial velocity distribution of 
the electron and its excursion time. The resonant enhancement is observed only for the long trajectory contribution of a particular 
harmonic when a window resonance in argon, which is off-resonant in the field-free case, is shifted into resonance due to a large 
dynamic Stark shift. 


The process of high-order harmonic generation (HHG) 
driven by a strong infrared (IR) laser field interacting with a 
rapidly ionizing medium is the main light source for the field 
of attosecond science The HHG process can be used to 
produce attosecond pulses because there is a natural, sub¬ 
cycle electron dynamics built into the physics of HHG^’^, 
which leads to a very broad plateau of emitted harmonics. 
This means that studying the HHG process itself in detail 
can, in principle, provide a deeper understanding of strong 
field electron dynamics at the attosecond time scale. Over 
the last decade experiments have indeed shown that the 
sub-cycle dynamics of HHG are encoded in the harmonic 
spectrum, though extracting them is complicated because of 
the highly non-linear nature of the process. 

Much of the promise in using HHG to better understand 
strong field physics at the sub-cycle level can be attributed 
to the effectiveness of the simple, semi-classical three-step 
model commonly used to describe the generation process 
In this model, an electron is first tunnel ionized and then ac¬ 
celerated by a strong laser field. If the electron is driven back 
to the vicinity of the ion by the oscillating strong field, the ac¬ 
cumulated energy may be emitted as a photon when the 
electron and ion recollide. The sequence of ionization and re¬ 
turn times leading to a specific harmonic frequency is loosely 
referred to as a trajectory because much of the physics can be 
understood by considering classical electron trajectories in a 
strong laser field, ignoring atomic effects after the ionization 
step and before the return. Depending on when during the 
laser cycle the ionization occurs, the electron will have differ¬ 


ent excursion and return times to the ion leading to different 
photon emission frequencies, resulting in a comb of odd har¬ 
monics of the laser field if the process is repeated over many 
laser cycles. Even in this simple model, however, there is not 
a one-to-one correspondence between the harmonic emission 
strengths and specific trajectories, because there are different 
trajectories leading to the same final energy. Trajectories that 
lead to the same photon energy interfere at the single atom 
level. The effect of this can either be studied or cir¬ 
cumvented, for example, through phase matching or spatial 
separation in the far field of the harmonics as is done in the 
present work. 

Studying trajectory resolved contributions to the HHG 
spectrum is an attractive proposition because different tra¬ 
jectories probe very different ionization conditions and have 
different excursion times. The most prominent contributions 
to the harmonic emission strengths come from the so-called 
short and long trajectories, which have excursion times of 
less than one laser cycle. Within a laser cycle the long tra¬ 
jectories are ionized close to the peak field strength and have 
an excursion time exceeding 0.65 laser cycles. The short tra¬ 
jectories are ionized at low field strengths and have shorter 
excursion times. Fortunately, the emission from these two 
trajectory classes can be separated experimentally in the far 
field enabling, for each harmonic frequency, comparison be¬ 
tween ionization at two different sub-cycle field strengths, 
followed by two different excursion times. This requires, 
ideally, that high accuracy measurements of both long and 
short trajectory contributions to each harmonic be made in 
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the same experimental setup 

Until recently most experimental efforts that make use of 
high harmonics have been concentrated on optimizing HHG 
from short trajectories, since their emission is well collimated 
and spectrally narrow. In addition, trains of attosecond pulses 
have been successfully created and measured by selecting 
the short trajectories’ contributionsThe emission from 
the long trajectories is more challenging to use because it 
is spectrally broader and more divergent, hence it is usu¬ 
ally removed by spatial filtering and/or the selection of spe¬ 
cific phase matching conditions in experiments as they oth¬ 
erwise can affect the temporal structure of the attosecond 
pulses In this paper we report on measurements made 
with very well-controlled, high repetition-rate laser pulses, 
which allow us to make trajectory resolved HHG measure¬ 
ments in argon gas while varying the ellipticity and the peak 
field strength of the driving laser pulses. The results allow 
us to elucidate new features in the sub-cycle ionization step 
that lead to long trajectories, that is, ionization at high field 
strengths followed by long excursion times. 

In this article, we present two methods of probing sub¬ 
cycle strong-field dynamics by comparing the trajectory re¬ 
solved emission of high harmonics and then studying the 
long trajectories in depth. In the first part of the article, a 
detailed experimental comparison of the ellipticity depen¬ 
dence as a function of harmonic order is presented, for both 
the short and the long trajectories. While harmonic gener¬ 
ation using elliptically polarized driving fields has been ex¬ 
tensively studied for the short trajectories both experimen¬ 
tally and theoretically the polarization dependence of 
the long trajectories have so far only been investigated theo¬ 
retically It follows from the simple three-step model 

that harmonic generation will be very sensitive to the elliptic¬ 
ity of the driving laser since the field acting on the electron 
while it is far from the ion can cause it to miss the recolli¬ 
sion. Since the long and the short trajectories have different 
excursion times, the impact of changing the ellipticity will 
be different for the two classes of trajectories. 

To explain the ellipticity dependence of the short trajec¬ 
tories it is sufficient to include wave packet spreading due 
to quantum diffusion, which we can model by including a 
distribution of momenta transverse to the instanteneous field 
vector at the moment of ionization. This distribution does not 
need to depend in detail on the moment of ionization, since 
the ionization field strength is low for short trajectories. In 
order to explain the ellipticity dependence of the long trajec¬ 
tories, however, this simple quantum diffusion model is not 
enough. Due to a larger variation in ionization field strength 
for the different long trajectories, a field-strength dependent 
momentum distribution has to be taken into account. We ex¬ 
pect that at higher field strengths a broader transverse mo¬ 
mentum distribution results from the lowering of the ioniza¬ 
tion barrier. We include this effect in our theoretical analysis 
of the long trajectory data via a simple extension of the three 
step model and find that it fits the our experimental data very 
well over the HHG plateau. 

In the second part of our study the sub-cycle sensitivity of 
trajectory resolved HHG measurements is used to study a re- 
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Figure 1: Typical harmonic spectra optimized to generate 
harmonics from both the short and the long trajectories. 

gion of the spectrum in which atomic resonances can alter the 
HHG signal. In particular, a window resonance in argon that 
is far from any harmonic of the laser frequency in the field- 
free case is shown to have a large effect on the long trajec¬ 
tory harmonic closest to it, but little or no effect on the short 
trajectory. We attribute this to the fact that the long trajec¬ 
tory component is dynamically Stark shifted into resonance 
by the laser field, which leads to a drastic enhancement of 
the emission from the long trajectory, but not the short where 
the field strength is much weaker and is not sufficient to shift 
the state into resonance. We measure this effect for a set of 
resonant harmonics over a range of driving field intensities. 

Experimental setup 

The experimental setup used for the experiment presented in 
this article is described in a recent publication^^ and is briefly 
outlined here. An YbiKGW based laser system (“Pharos”, 
Light Conversion Ltd.) was used to deliver 170 fs, pulses 
with a central wavelength of 1030 nm. The laser system has 
a variable repetition rate between 1 and 600 kHz, but all the 
presented data were recorded at a repetition rate of 20 kHz. 
The pulses were focused tightly into a continuous argon gas 
jet, with a 90 \im orifice, using a 100 mm focal length achro¬ 
matic lens. Directly after the interaction region, a differential 
pumping hole with an inner diameter of 0.5 mm was placed 
to minimize the background gas in the detection chamber. 
The differential pump hole allowed for a pressure difference 
of the background gas between the generation and detec¬ 
tion chambers of 4-5 orders of magnitude. The HHG spec¬ 
trum was measured using a home-built imaging spectrometer 
based on a variable-line-spacing grating and a microchannel- 
plate with an attached phosphor screen and a camera with a 
resolution of 2456x2058 pixels and a dynamic range of 14 
bits. The grating diffracts and refocuses the XUV in the hori¬ 
zontal direction while the vertical direction is left unaffected. 
Therefore the vertical direction provides the divergence of 
the XUV light while the horizontal direction shows the spec¬ 
trum. 

Ellipticity measurement 

Figure 1 shows a typical harmonic spectrum when a linearly 
polarized driving laser is used and the gas jet is placed in 
the focal plane of the generating beam. The experimental 
parameters (pulse energy, gas density, spot size, etc.) were 
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optimized to generate harmonics from both the short and the 
long trajectories. 

The Gaussian transverse and temporal intensity profile of 
the driving laser, in combination with the fact that the dipole 
phase of the long trajectories has a stronger intensity depen¬ 
dence than the short trajectories, result in a larger wavefront 
curvature and more divergent light generated by the long tra¬ 
jectories this also explains the spatial-spectral rings ob¬ 
served in the far field. We therefore attribute the inner part 
of the harmonic spectrum to be dominated by the short tra¬ 
jectories, while the spatial-spectral rings are attributed to in¬ 
terference between long trajectories of different emitters. As 
the trajectory dependent dipole phase is strongest for the low 
orders, the interference rings are mainly seen for the low end 
of the plateau region. This spatial separation was exploited 
to study the contributions from the long trajectories only. 

A quarter-wave plate was used to introduce ellipticity, de¬ 
fined as the ratio between the minor and major axis compo¬ 
nents of the driving laser field. Figure 2 (a) shows an en¬ 
larged view of the spatial-spectral profile of harmonics 19- 
23 of Fig. 1. Figure 2 (b) shows a measurement of the spa¬ 
tially and spectrally integrated strength of harmonic 23 (H23) 
as a function of ellipticity of the driving laser relative to the 
strength at linear polarization. The integrated signal clearly 
follows a Gaussian distribution with respect to ellipticity as 
previously observed 

We define the threshold ellipticity, as the amount of 
ellipticity required for the harmonic signal to drop by a fac¬ 
tor of two compared with linear polarization. Our very high 
signal-to-noise ratio allows us to analyze the ellipticity de¬ 
pendence of each pixel rather than the spatially and spec¬ 
trally integrated signal. Figure 2 (c) presents the strength 
of three different pixels within H23 as a function of ellip¬ 
ticity. The strength of each pixel is fitted with a Gaussian 
profile to extract the corresponding threshold ellipticity of 
each pixel, which are used to create a two-dimensional map 
of the threshold ellipticity as a function of energy and di¬ 
vergence angle. The full threshold ellipticity maps for the 
conditions of Fig. 1 can be found in the Methods section. 
Focusing on H23 in Fig. 2 (d), we observe three different re¬ 
gions of threshold ellipticity; an inner region with a threshold 
ellipticity around 0.16, and two outer regions with threshold 
ellipticities of around 0.09 and 0.1 respectively. 

Figure 3 presents the average threshold ellipticity for both 
the long and the short trajectories as a function of harmonic 
order for the conditions of Fig. 1. For the short trajectories, 
we observe a behavior similar to previous measurements ^ 0 ^ 22,26 
i.e., the threshold ellipticity decreases slowly with increasing 
harmonic order. A similar trend is also observed for the long 
trajectories, in contrast to what would be expected if only the 
excursion time of these trajectories is considered. 


Ellipticity theory 

In a semi-classical model where the propagation step is cal¬ 
culated classically and the electron only has a velocity par¬ 
allel to the electric field, only linearly polarized light would 
produce high-order harmonics since any ellipticity will pre- 
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Figure 2: (a) Enlarged view of the harmonic spectrum at lin¬ 
ear polarization for H19-H23 of Fig. 1. (b) Spatially and 
spectrally integrated signal of H23 as a function of elliptic¬ 
ity. (c) Measurement of the ellipticity dependence of three 
different spatial-spectral parts of H23 as indicated by the 
three arrows. The solid lines represent Gaussian fits to the 
experimental data, (d) Pixel-by-pixel threshold ellipticity of 
the spatial-spectral region of part (a). 



Figure 3: The red (blue) line show the measured threshold 
ellipticity as a function of harmonic order for the long (short) 
trajectories when the laser is focused in the middle of the gas 
jet. The error bars indicate the standard deviations of the 
corresponding threshold ellipticity. 


vent the electron from returning to its original position. The 
fact that high-order harmonics are observed even for ellipti- 
cally polarized light is usually attributed to quantum diffu¬ 
sion; the electron wave packet spreads out as it is accelerated 
in the laser field. The wave packet spread allows for an over¬ 
lap between the electron and the parent ion, even when the 
electron is transversely displaced due to the elliptically po¬ 
larized laser field. 

Quantum diffusion can be seen as resulting from an ini¬ 
tial distribution of velocities of the electron - the more con¬ 
fined the electron is in one direction, the more it will spread. 
In particular, a spatial confinement in the direction perpen¬ 
dicular to the laser field, will lead to a transverse velocity 
distribution which is necessary for HHG. A rough estimate 
of the confinement is given by the size of the groundstate. 
Using this estimate, a trajectory spending longer time in the 
continuum will diffuse more which results in a lower HHG 
yield. 

For the short trajectories, the above estimate of the quan¬ 
tum diffusion, which is independent of the ionization time, is 
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sufficient to explain the increase in sensitivity, as a function 
of harmonic order. For this set of trajectories the highest en¬ 
ergy photons are produced by electrons with the longest ex¬ 
cursion time. As the transverse displacement of the electron 
at the point of recombination increases with the excursion 
time, the trajectories leading to the higher harmonics are dis¬ 
placed more than those leading to the low orders. Therefore 
the overlap between the ion and the electron at the recombi¬ 
nation time decreases with harmonic order. 

For the long trajectories this effect leads to the opposite 
result, as the kinetic energy of the returning electrons de¬ 
creases with increasing excursion time. To understand the 
experimentally observed ellipticity dependence of this set of 
trajectories, we apply a model that also takes the sub-cycle 
variation of the initial electron velocity distributions into ac¬ 
count, as well as the change in excursion time for the differ¬ 
ent trajectories as the ellipticity is varied. This effect plays 
a major role in the initial velocity distribution as the long 
trajectory electrons are ionized closer to the peak of laser 
field, where the atomic potential is more distorted in the di¬ 
rection of the laser field, and thus the electron wave packet 
is more perpendicularly confined at the time of ionization 
The perpendicular confinement of the electron at the ioniza¬ 
tion time leads to a large uncertainty in the perpendicular 
velocity distribution. 

Our method is similar to ^4-26^ definition of thresh¬ 

old ellipticity is not the same in the different studies. The 
procedure is as follows: First, we calculate the return en¬ 
ergy of the electrons for the two first sets of trajectories as 
a function of both ionization time and ellipticity. The po¬ 
sition of an electron released at time ti in an elliptical field 
F{t)=F/VlT^ [sin((or);ecos((or)] is found by integrat¬ 
ing the Newtonian equations of motion twice: 


. s _ F r sin(( 0 r) - sin{coti) - co{t - u) cos{coti) 
~ m^Vl + e^ e[cos((or) - cos{coti) + co{t - ti) sin((or/)] 

F{t — ti)\i + r/, 

( 1 ) 

where and are the initial velocity and position, respec¬ 
tively; F is the field amplitude, (0 the frequency of the fun¬ 
damental field, e G [—1, -hi] is the ellipticity, with 0 meaning 
linear polarization along the x axis. Atomic units are used. 
We assume that = r(r^) = 0, where ty is the moment of 
return. Finding this time requires solving the transcendental 
equation numerically. For elliptical polarization, the drift ac¬ 
quired by the electron can be countered by an initial velocity 
V; that is transverse to the driving field at the time of ioniza¬ 
tion (this is analogous to quantum diffusion of the electron 
wavepacket as it is accelerated in the laser field). Thus, we 
solve (1) for ty and v/ for each U G [0.257,0.57], 7 being the 
period, and each £ G [0,1]. We assume that \i = V|| + v^, 
where the two components are parallel and perpendicular to 
the driving field at the time of ionization. Furthermore, we 
assume that V|| =0, such that all uncertainty is in the initial 
transverse momentum, = meV^, {mq = 1 in atomic units). 

The kinetic energy at the time of return is given by W[ = 
p^{ty)j2\ this gives the map of energies seen in Fig. 4. As 



Figure 4: Map of return energies, as a function of elliptic¬ 
ity £ and ionization time ti. /p is the ionization potential of 
the ground state. Plotted are also isoenergetic curves cor¬ 
responding to the harmonics of the fundamental field, and 
the cut-off energy (shown by the solid black line), which de¬ 
creases for increasing ellipticity. Trajectories, which are ion¬ 
ized earlier than the cut-off energy, correspond to the long 
trajectories. It is easy to see that a specific time of ionization 
does not correspond to a certain recombination energy. 


can be seen in Fig. 4, the cut-off position is shifted when the 
ellipticity is increased {i.e. the highest energy photons can 
only be produced from linearly polarized light) and the initial 
timing leading to a specific harmonic order is also changed. 
This trend is even more clear when lineouts at different ellip- 
ticities are presented as in Fig. 5. 

The next step to estimate the harmonic yield is to calcu¬ 
late the combined probability of ionizing at time ti and having 
the required initial velocity for the electron to return. This 
is possible since the correspondence between a certain har¬ 
monic and its ionization time for different ellipticities is al¬ 
ready calculated. Since the required transverse momenta for 
the long trajectories to return are quite large for high elliptic¬ 
ities, we use the full expression for the transverse momentum 
distribution found in reference 33, 




2F{ti) 


( 2 ) 


/p is the ionization potential of the ground state. The tunnel¬ 
ing rate is taken from ADK theory This combined prob¬ 
ability, which is the product of the separate probabilities de¬ 
scribed above, is displayed using a colour scale in Fig. 6, as 
a function of ellipticity and ionization time. The isoenergetic 
curves from Fig. 4 are also included in the figure. To cal¬ 
culate the yield of a given harmonic order as a function of 
ellipticity, one extracts the probability along the correspond¬ 
ing isoenergetic curve. 

Finally, in the last step of our model, the calculated yield 
as a function of ellipticity is fitted with Gaussian functions 
for each harmonic to obtain the threshold ellipticities in a 
similar manner to the experimental data. The result of the 
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Figure 5: (a) Lineouts of Fig. 4 for 11 equidistant ellipticities 
from linear to circular polarization. From this plot it is obvi¬ 
ous that the cut-off decreases with increasing ellipticity and 
that it occurs for earlier ionization times. It is also clear that 
the initial timing necessary to produce a specific harmonic 
changes with the ellipticity. (b) The same as (a), but plotted 
as a function of excursion time te=tr — U instead. 


model for the long trajectories is presented in Fig. 7 (a) to¬ 
gether with the experimental data. 

Our extended model compares very well with the exper¬ 
iment presented in this work; in particular the decrease in 
threshold ellipticity with increasing harmonic order is ex¬ 
plained. This is opposite to what would be expected from 
sub-cycle field-independent quantum diffusion. It is also op¬ 
posite to the analytical expression presented in reference 25 
which is included for comparison in Fig. 7 (a) as a dotted 
line. 

From Fig. 7 (a) it is clear that some of the long trajectory 
harmonics (17, 21 and 23) have a lower threshold elliptic¬ 
ity than what is predicted by the model. We attribute this 
to the presence of atomic resonances in the vicinity of the 
corresponding energies, which clearly cannot be captured by 
the model we are using. In what follows, we demonstrate 
that these resonances can be dynamically Stark shifted by 
the sub-cycle field strength and will therefore infiuence the 
harmonic generation differently for the short and the long 
trajectories. The spatial separation of the short and long tra¬ 
jectories leading to the same energy, enable us to directly 
compare the infiuence of the sub-cycle field strength. For 
a given harmonic order, the short trajectory is initiated at a 
field strength which is insufficient to shift the state into reso¬ 
nance, thereby precluding the enhancement observed for the 
long trajectory initiated at a higher sub-cycle field strength. 


Resonant HHG 

Resonant HHG in argon at a photon energy corresponding to 
H17 in our experiment has previously been observed In 
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Figure 6: Map representing the combined probability of ion¬ 
izing at time ti and having the initial transverse velocity re¬ 
quired for return as a function of ellipticity and ionization 
time. The isocurves are the same as in Fig. 4, represent¬ 
ing constant return energy. Following an isocurve gives the 
probability of generating a certain harmonic, as a function of 
ellipticity. 


our study, we also see an effect for H21 and H23 [Fig. 7 (b)]. 
We will focus the discussion on H21, where the effect is most 
pronounced. At linear polarization a strong enhancement of 
H21 is clearly observed, while this enhancement is gone for 
an ellipticity of 0.2 as can be seen in Fig. 7 (b). The change in 
ellipticity leads to a variation in the intensity, since the pulse 
energy is kept constant. This means that the observed effect 
can be due to either the intensity or the ellipticity. In order 
to disentangle the two, an intensity scan was performed for 
linearly polarized light. Figure 8 shows the experimentally 
measured intensities of the long (a) and short (b) trajectories 
for H15-H37 as a function of IR pulse energy. As clearly ob¬ 
served, the yield of H21 rises more rapidly for the long trajec¬ 
tories once the IR pulse energy exceeds 0.14 mJ [Fig. 8 (a)], 
while the emission from the short trajectories are left unaf¬ 
fected [Fig. 8 (b)]. Enhancement of the long trajectories can 
also be clearly seen for H17 and H23, albeit, at slightly lower 
pulse energies. Harmonic 17 qualitatively follows the trend 
predicted for the single atom response given in reference 37, 
thus H17 is not further discussed in the present work. Har¬ 
monic 23 reaches a maximal strength at 0.15 mJ whereafter 
a slow decrease with respect to increased pulse energy is ob¬ 
served. 

We interpret the behavior of H21 (but also H17 & H23) 
to be the result of HHG in the presence of an atomic res¬ 
onance. Resonant HHG may increase the harmonic yield 
through a number of different mechanisms . For H21 
with photon energy of 25.2eV, the closest resonance is the 
3s^3p^ -> 3s^3p^4p^ transition (26.6 eV), which is a window 
resonance . Our interpretation requires that the 3s^ 3p^4p^ 
state, which is lowest state in the 3s-^np closed channels, is 
red-shifted by approximately 1.4 eV [see inset in Fig. 8 (a)]. 
This is feasible as the dynamical Stark shift of the 3s^3p^4p^ 
state should be dominated by the interaction with the 3s-^np 
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Figure 7: (a) The dashed lines show the numerically calcu¬ 
lated threshold ellipticities for the long trajectories in the 
plateau region according to our model. Correspondingly, 
the dotted lines show the threshold ellipticity given by ref¬ 
erence^^. For comparison, the experimental data for the long 
trajectories are shown. The simulations were done for a laser 
intensity of 8.5 • 10^^ W/cm^ and a wavelength of lOSOnm. 
(b) Integrated harmonic spectra for in-focus generation with 
various ellipticities. The color scale corresponds to differ¬ 
ent ellipticities from black (e = 0) to light gray (e = 0.25) 
in steps of 0.05. The fillings between the different harmon¬ 
ics allow us to better visualize the differences. The various 
spectra have been normalized to the level of HI5. The error 
bars used in (a) show the standard deviations of the threshold 
ellipticities of the corresponding harmonic. 


closed channel, rather than through coupling with the 3p^ 
Z(s,d) open channelsand the Stark shifts on the order of 
the ponderomotive energy are well known In addition 
to the Stark effect, the IR intensity also causes a blueshift of 
the IR energy, and thereby the XUV photon energies. How¬ 
ever, it was confirmed from the intensity scan that this effect 
is too small to explain the results, as the central frequencies 
of the harmonics did not change. Since we only observe the 
enhancement for the long trajectory, we interpret this as an 
effect of the comparably higher field strength for this trajec¬ 
tory, at the time of ionization. 

Reshaping of the argon HHG spectrum by this particular 
window resonance has previously been observed however, 
the interpretation is fundamentally different in the work pre¬ 
sented here. In reference 48, a few-cycle pulse was used 
to generate broadband harmonics by HHG in a gas jet. The 
backing pressure for the continuous gas jet was then increased 
significantly so that all other wavelengths than exactly the 
resonant wavelength were suppressed by re-absorption in the 
generating gas. This led to a narrowing of H17 (of 800 nm) 
from a width of roughly 1.5 eV to a width comparable to the 
field-free width of the window resonance This leads us to 
believe that the effect observed in reference 48 happens over 
a large volume, where the IR intensity is weak. The results 
presented in this article, however, is clearly an effect that take 
place at high laser intensity, where the dynamical Stark effect 
is strong. 

The behavior of H23 with respect pulse energy of the 
driving laser can be understood as an effect of over-shifting 
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Figure 8: Experimentally measured long (a) and short (b) 
trajectory intensities of H15-H37 as a function of pulse en¬ 
ergy of the fundamental field. H17 is shown in gray, H21 in 
black, and H23 in brown, while the remaining harmonics are 
shown in a rainbow color scale from dark blue (HI5) to dark 
red (H37) with increasing order. The dashed line indicates 
roughly the conditions of Fig. 1. The inset of (a) schemat¬ 
ically shows how the Rydberg series for the 3 s electrons in 
combination with the continuum for the 3p electrons creates 
a series of window resonances. Absorption of an XUV pho¬ 
ton can create a coherent superposition of the two valence 
electrons which interfere and affects the absorption cross- 
section (Tabs • In argon this amounts to a reduction of the ab¬ 
sorption cross-section^^. The vertical arrows within the inset 
indicate the expected direction of the light induced energy- 
shift of the respective 3s3p6np states. 
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of the atomic resonances causing the enhancement. As indi¬ 
cated in the inset of Fig. 8 (a) the field-free detuning of H23 
is less than the detuning of H21 with respect to both the tran¬ 
sition energies into the 3s^3p^4p^ and the 3s^3p^5p^ states, 
so the minimal required pulse energy for enhancement of this 
harmonic is lower. Nevertheless, the maximal enhancement 
factor is largest for H21, due to the strong dipole coupling 
with the red-shifted 3s^3p^4p^ state. 

As the pulse energy is increased beyond the optimum en¬ 
ergy for resonant generation of H23 the enhancement of this 
harmonic starts to vanish. The slow decrease likely origi¬ 
nates from the long pulse duration of the driving laser, which 
means that a number of cycles will have the optimum energy 
shift. A similar effect is expected to occur for H21 at higher 
pulse energies, however, due to limitations of the laser sys¬ 
tem this was not seen in the present work. 

Apart from the major effects on H17, H21 and H23 ob¬ 
served both in the ellipticity and the intensity measurements 
a minor amount of enhancement of HI9, H25 and H27 can 
be observed for the long trajectories once the IR pulse energy 
exceeds 0.15 mJ [Fig. 8 (a)]. The field-free detuning from the 
3s3p6np manifold of resonances are larger for these harmon¬ 
ics, so any enhancement effect on these harmonics is both 
expected to be less, and to occur at higher pulse energies in 
full agreement with the observation. 

In conclusion, we have experimentally investigated the 
ellipticity and intensity dependencies of HHG from the long 
and the short trajectories. This type of measurements en¬ 
ables us to probe the infiuence of the sub-cycle field strength 
on HHG process. We have shown that the well-established 
semi-classical model has to be extended by taking the instan¬ 
taneous field strength into account, to also describe the gen¬ 
eral behavior of the long trajectories. We have demonstrated 
how off-resonant states embedded in the continuum can en¬ 
hance long trajectory harmonics by being shifted into reso¬ 
nance by the strong driving laser, different amounts for dif¬ 
ferent trajectories due to the sub-cycle nature of the genera¬ 
tion process. When the driving laser field is strong enough to 
cause an enhancement at linear polarization, these harmon¬ 
ics show a stronger ellipticity dependence as the dynamical 
Stark shift depends on the polarization. 

This study highlights the importance of systematical stud¬ 
ies of the generation process under various conditions. Fur¬ 
thermore, the extension of the knowledge of the harmonic 
generation process to the long trajectories will be beneficial 
for high-order harmonic spectroscopy studies. 
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Methods 


Evaluation of experimental data 

In this section we present the details of the analysis method used for 
the experimental data. 

It is well-known from the strong field approximation that for har¬ 
monics in the plateau region, there are several electronic trajecto¬ 
ries, which may contribute to the generation process. Emission from 
these different trajectories interferes and shapes the far field spatial 
spectral profile. The phase of these trajectories can be approximated 
with a phase proportional to the intensity I{x,y,zj) such that 

0^> = a^h{x,y,z,t), (3) 

where (j)q^^ is the trajectory dependent dipole phase and is the 
proportionality constant. The first two sets of electron trajectories 
are usually referred to as the short and long trajectories. It well- 
established that in the plateau region the proportionality constants 
are much larger for the long trajectories than for the short trajecto- 


The short trajectories can be isolated in the generation process by 
placing the gas jet behind the focal plane and adjusting gas pres¬ 
sure and pulse energy accordingly A spectrum optimized for 
this is shown in Fig. 9 (a). Figure 9 (b) shows the corresponding 
threshold ellipticity map, which is extracted in a similar manner as 
in the main article. When the gas jet instead is placed at the focus of 
the laser both sets of trajectories can efficiently be phase-matched 
by adjusting the other experimental parameters accordingly. Fig¬ 
ure 9 (c) shows a spectrum optimized to generate with both sets of 
trajectories, while Fig. 9 (d) is the corresponding threshold elliptic¬ 
ity map. As a consequence of the larger dipole phase of the long 
trajectories, the light produced by the these trajectories are more di¬ 
vergent. This effect was used to spatially separate the contributions 
from only the long trajectories in Fig 9 (c,d). 

Figure 10 shows enlarged views of H19 for out-of-focus genera¬ 
tion [(a)] and in-focus generation [(b)]. In the out-of-focus case the 
harmonic exhibits a homogeneous spatial-spectral dependence with 
respect to ellipticity, this is not observed for the in-focus case, where 
several regions of ellipticity dependence are clearly observed. 

The homogeneity of the ellipticity dependence for out-of-focus gen¬ 
eration [Fig. 10 (a)] reveal that in order to study ellipticity depen¬ 
dence of the short trajectories out-of-focus an imaging spectrometer 
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Figure 10: (a) An enlarged view of the ellipticity map of 
H19 for out-of-focus generation [Fig. 9 (b)]. The lower 
panel shows a histogram of the threshold ellipticities. Each 
pixel is weighted with the corresponding pixel strength at 
linear polarization, (b) An enlarged view of the ellipticity 
map of H19 for in-focus generation [Fig. 9 (d)]. The lower 
panels show histograms of the threshold ellipticities within 
the white boxes (left panel) and the black box (right panel). 
The histograms are weighted with the corresponding pixel 
strength at linear polarization. 


Figure 9: (a) Harmonic spectrum at linear polarization when 
the gas jet is placed behind the focal plane of the laser, (b) 
Threshold ellipticity map as a function of energy and diver¬ 
gence angle for out-of-focus generation, (c) Harmonic spec¬ 
trum at linear polarization when the gas jet is placed in the 
focal plane of the laser, (d) Threshold ellipticity map as a 
function of energy and divergence angle for in-focus genera¬ 
tion. 


Detection efficiency 

We measured the ratio in detection efficiency between horizontal 
and vertical polarization to be 1.38. 3D-TDSE calculations^^ show 
that the plateau harmonics exhibit a smaller ellipticity than the driv¬ 
ing infrared laser. We therefore estimate that the upper limit of the 
non-fixed major axis impact on the measurement to be given by the 
following expression: 


is not needed and a spatial-spectral integration with respect to har¬ 
monic order would be sufficient. This is clearly not the case of 
in-focus generation [Fig. 10 (b)]. 

Figure 11 shows the threshold ellipticity as function of harmonic or¬ 
der for the conditions of Fig. 9. Figure 11 shows normalized thresh¬ 
old ellipticity histograms as function of harmonic order. In the in¬ 
focus generation case both the on-axis emission and off-axis emis¬ 
sion are shown, in (b) and (c) respectively, while for out-of-focus 
case only the on-axis emission is shown in (a). We note that for har¬ 
monics close to the cut-off indications of the long trajectory appears 
also for the out-of-focus case. In order to extract expectation values 
and standard deviations for the various threshold ellipticities as a 
function of harmonic order and trajectories the experimental data is 
smoothened using the Kernel density estimation methodAfter 
smoothing the data was fitted with two Gaussian distributions for 
the long trajectories, while the short trajectories where fitted with 
a single Gaussian distribution. The expectation value of the fitted 
Gaussian distributions are plotted as solid lines in Fig. 11, while the 
uncertainty bars show the corresponding standard deviations of the 
fits. 


fdet DaDfy , (4) 

where Ea and Ejy are the major and minor axis component of the 
infrared electrical fields, and Da and D^y are the respective detection 
efficiencies. Using this expression together with the standard Jones 
matrix calculus for polarization of the infrared light we estimated 
the upper limit on the determination of the threshold ellipticity to be 
less than the presented standard deviations. The presented data was 
performed around the linear polarization direction with the highest 
detection efficiency. Therefore the threshold ellipticity might be 
systematically underestimated slightly. 
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